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We consider a cavity optomechanical cooling configuration consisting of a mechanical resonator 
(denoted as resonator 6) and an electromagnetic resonator (denoted as resonator a), which are 
coupled in such a way that the effective resonance frequency of resonator a depends linearly on the 
displacement of resonator b. We study whether back-reaction effects in such a configuration can be 
efficiently employed for suppression of decoherence. To that end, we consider the case where the 
mechanical resonator is prepared in a superposition of two coherent states and evaluate the rate of 
decoherence. We find that no significant suppression of decoherence is achievable when resonator a is 
assumed to have a linear response. On the other hand, when resonator a exhibits Kerr nonlinearity 
and/or nonlinear damping the decoherence rate can be made much smaller than the equilibrium 
value provided that the parameters that characterize these nonlinearities can be tuned close to some 
specified optimum values. 

PACS numbers: 



I. INTRODUCTION 

The quest for quantum effects in nanomechanical de- 
vices has motivated an intense research effort in recent 
years [l|-|3| . Experimental demonstration of quantum su- 
perposition in a nanomechanical resonator may provide 
an important insight into the problem of quantum to clas- 
sical transition |4l4lOl|. However, in many cases the life- 
time of such superposition states is too short for experi- 
mental observation since the coupling between a nanome- 
chanical resonator and its environment typically results 
in rapid decoherence [ll|, [l^] ■ As a case study, consider 
a superposition of two coherent states \ai) and |a2) of 
a mechanical resonator having an angular resonance fre- 
quency iOb and damping rate 7f,. The decoherence rate of 
such a superposition state is given in the high tempera- 
ture limit A:bT > fiwfe by [H-al) 



— =47h (5a - — 



(1) 



where 5a — a2 ~ ai. 

While Eq. ([T]) was derived by assuming linear response, 
it is well known that nonlinear response can be exploited 
for reduction of thermal fluctuations. One example is 
the technique of noise squeezing that can be employed 
for reducing thermal fluctuations in one of the quadra- 
tures of a mechanical resonator [13, [3 ■ Another exam- 
ple, which is the focus of this chapter, is the techiiique of 
optomechanical cavity cooling. This technique |19l426l |. 
which was first proposed as a way to enhance the detec- 
tion sensitivity of gravity waves [23, ^^ , can be employed 
for significantly reducing the energy fluctuations of a me- 
chanical resonator well below the equilibrium value [29l - 
|43 |. Cooling is achieved by coupling the mechanical res- 
onator (denoted as resonator b) to an electromagnetic 
resonator (denoted as resonator a) in such a way that 
the effective resonance frequency of resonator a becomes 
linearly dependent on the displacement of resonator b. 
When the parameters of the system are optimally cho- 



sen the fluctuations of resonator b around steady state 
can be significantly reduced well below the equilibrium 
value by externally driving resonator a with a monochro- 
matic pump tone. In this region back-reaction due to 
the retarded response of the driven resonator a to fluc- 
tuations of resonator b acts as a negative feedback, pro- 
viding thus additional damping which results in effective 
cooling down of resonator b. The success of these exper- 
iments raises the question whether similar back-reaction 
effects can also be efficiently employed for suppression of 
decoherence below the equilibrium value. 

Here we study this problem by generalizing Eq. ^ 
for the case where cavity cooling is applied. Nonlinearity 
in resonator a is taken into account to lowest nonvanish- 
ing order. The equations of motion of the system are 
obtained u sing the Gardiner and CoUett input-output 
theory |43l |44| . By linearizing these equations we de- 
rive the susceptibility matrixes of the system, which al- 
low calculating the response of both resonators to input 
noise. This, in turn, allows evaluating both, the spectral 
density of fluctuations and the decoherence rate l/r,^ of 
resonator b. In both cases we examine the cooling effi- 
ciency by defining an appropriate effective temperature 
and by calculating it for an optimum choice of the sys- 
tem's parameters. We find that only modest suppression 
of decoherence is possible using cavity cooling unless the 
system is driven into the region of nonlinear oscillations. 



II. THE MODEL 

The model consists of two resonators, labeled as a and 
b respectively, which are coupled to each other by a term 

hflNa (Ab + Al) in the Hamiltonian. Here Aa, A^^ and 

Na = AlAa {Ab, AI and Nb = AlAb) are respectively 
annihilation, creation and number operators of resonator 
a (6). The first resonator is coupled to 3 semi- infinite 
transmission lines. The first, denoted as a 1, is a feedline, 
which is linearly coupled to resonator a with a coupling 



constant Tai, and which is employed to deliver the input 
and output signals; the second, denoted as a2, is linearly 
coupled to resonator a with a coupling constant Ta2, and 
it is used to model linear dissipation, whereas the third 
one, denoted as aS, is nonlinearly coupled to resonator a 
with a coupling constant Tas, and is employed to model 
nonlinear dissipation. Linear dissipation of resonator b 
is modeled using semi-infinite transmission line, which is 
denoted as b and which is linearly coupled to resonator 
b with a coupling constant Tf,. Kerr- like nonlinearity of 
the driven resonator a is taken into account to lowest 
order by including the term {h/2) KaA'lAj^AaAa in the 
Hamiltonian of the system, which is given by 



H - huJaNa + -KaAlAlAaAa + hWbNb 



hnNa A. 



duj a\^ (w) aai (w) ui 



AljJ 



T,i4a,iH-fT:ialiHA, 



dw ajj2 (w) aa2 (w) w 



AljJ 



Ta2Alaa2 (uj) + T*^a\^ (w) Aa 



Alo a\^ (w) aa3 (w) w 



Auj 



T,344«a3 (^) + T^galg {UJ) AaAa 



Aio a\ {io) Of, {ijj)uj 



Auj 



TbAlat (w) -t- ^al {uj) At 



(2) 



one has 
dAa 



-^ = -iuJaAa - iKaNaAa - inAa (^fc + AI 

-i Alu Taiaai {i^)-i / Auj Ta2aa2 (w) 



-2i Auj TasAluas (w) 



and 



dAb 

At 



-iujbAb — iflNa — i Auj Tbab (w) 
Using the bath modes commutation relations 

Gal (uj) , aj^i {uj') ^ S{uj -uj') , 



Iqi {uj) ,aai {uj')] = 



one obtains 



dagi {uj) 
At 



-iujUal (w) - JT^i^a 



(8) 
(9) 

(10) 
(11) 

(12) 



Using initial condition Uai {uj, tg) one finds by integration 
that 

aai(cj,t)=a„i(cj,to)e*"(*''-*) 

-iT*^ / di' Ag(i')e'"(*'"*) . 

Jto 

(13) 

Next we integrate Eq. (fT3|) over uj. The coupling co- 
efficient Tqi, which is assumed to be uj independent, is 
expressed as 



T -, — . / 7al i0„i 



(14) 



A. Equations of Motion 

The Heisenberg equations of motion are generated ac- 
cording to 



,n^ = p^n] , 



(3) 



where O is an operator. Using the commutation relations 



AbAl 



[Aa,Al] 

[Aa,Na\ = Aa, 
\Ab, A^fc] = Ab , 



1 , 



(4) 

(5) 
(6) 

(7) 



where 7ai is positive and 0ai is real. Using the following 
relations 



da;e'"(*'-*) ^2Ti6{t-t') , 



where sgn(x) is the sign function 

, , f +1 if X > 
sgn(.T)= ^ .^^^Q_ 



one finds that 
1 



d^ aai {uj,t) = <1 {t) - J j2|ie~'*»iA„ {t) 




(15) 



At'S{t-t')f{t')^-sgn{t-to)f{t) , (16) 



(17) 



(18) 



where 



<i it) 



1^ 



dc^aai(^,io)e'"^*°'*^ • 



(19) 



Using similar definitions the above results are general- 
ized for the other semi-infinite transmission lines that are 
linearly coupled (labeled as a2 and 6). For the transmis- 
sion line a3, which is nonlinearly coupled, the coupling 
coefficient Ta3, which is also assumed to be oj indepen- 
dent, is expressed as 



and the following holds 
1 



/ 7a3 i0^ 

27r 



(20) 



^ , dw aa3 (^, t) = <"3 (t)-jX2^e-*^- A, {£) Aa (i) 
27r J ^ 

(21) 
Substituting these results into Eqs. ([5]) and © yields 

dA, 



dt 



: - {iuJa + 7a + (^-ft^a + ^az) ^a] ^a 

2*V7^e'^-AtaS (t) , 



(22) 



and 



di 



: - (iwfc + 7b) Afc - il77Va 



where 



7a = 7al + 7a2 • 



B. Rotating Frame 



(23) 
(24) 



^-«-«' 



3 

(30) 



where 

"a == "a ( L/Q, C^, Ab, A(j 

= |i [Aa + O (a + 4)] + 7a + (iXa + 7a3) A^a } Ca 



-fiv/2^e*^°'&P- 



Aq = Wa - Wp , 



i^a = -iv/2^e**-Ci - iy2^e^^-C2 
- 2zV7^e^(^-+-p*)clC3 , 

= (iwfc + 76) Ab + i^Na 



(31) 
(32) 

(33) 



and 



Fb = -zv/2^e^*'4"(i) . 

III. LINEARIZATION 

Expressing the solution as 

^a = i>Q + Ca , 
Ab ^ Bb + Cb , 



(34) 
(35) 



(36a) 
(36b) 



where both Ba and Bb are complex numbers, and con- 
sidering both Ca and Cb as small one has to lowest order 

Bq (Ca,Ca,Cb,C^ ) = 9q {Ba,Ba,Bb,Bl) 



WiCa + W2CI + WsCb + WiC 



h ' 



(37a) 



Consider the case where a coherent tone at angular 
frequency Wp and a constant complex amplitude 6p is 
injected into the feedline. The operators of the driven 
resonator and its thermal baths are expressed in a frame 
rotating at frequency Wp as 



in _ in -iupt 



(25) 
(26) 
(27) 
(28) 



Using this notation Eqs. ([^ and (^5)) can be rewritten 
as 



dCa 

di 



Oa^Fa 



(29) 



©b (Ca,C],Cb,Cb — 9b {Ba,B*a,Bb,Bl) 



where 



+ W^Ca + WqcI + T^TCb + W^cl , 

(38a) 

M^l = iAf + 7a + 2 (iifa + 7a3) |Sa|' , (39a) 

T^2 = (ii^a+7a3)S2, (39b) 

W:i = Wi = inBa , (39c) 

T4^5 = if^S* , (39d) 

W^6 = i^Ba , (39e) 

W^7 = ioJb + 76 , (39f ) 

Ws = , (39g) 



and where 



Af^Aa + n{B, + B;) 



A. Mean Field Solution 



(40) 



Mean field solutions are found by solving 
©b {Ba,B*^,B]y,Bl) = , 



(41a) 
(41b) 



that is 



and 



iAf+-fa + {lKa+la:i)\Ba\'\Ba 



{iuJb + -fb)Bb + in\Bar =Q. 



(42) 



(43) 



Extracting Bi, from Eq. (|43p and substituting it in Eq. 
(|42l) yields 

[lAa +la+ {tKf + 7,3) \Ba\^}Ba 



+ i^2^e''^'^'bp = 0, 



where if° , which is given by 






(44) 



(45) 



is the effective Kerr constant. Taking the module squared 
of Eq. dUl) leads to 



{Aa+KfEa) +{la+lazEar Ea=2jalK 



where 



(46) 



(47) 



Finding Ea by solving Eq. (H51) allows calculating Ba 
according to Eq. (|i^ and Bf, according to Eq. (P5|. 



B. Onset of Bistability Point 



In general, for any fixed value of the driving amplitude 
bp, Eq. (|44)) can be expressed as a relation between Ea 
and Aa- When bp is sufficiently large the response of the 
system becomes bistable, that is Ea becomes a multi- 
valued function of Aq in some range near the resonance 



frequency. The onset of bistability point is defined as the 
point for which 



9Aa 

dEa 

d{Eaf 



= 



= 



(48) 
(49) 



Such a point occurs only if the nonlinear damping is suf- 
ficiently small [43| , namely, only when the following con- 
dition holds 



l^:*"! > V37a: 



3 • 



(50) 



At the onset of bistability point the drive frequency and 
amplitude are given by 



(Aa), 



K'= 



-la 



\Kf\ 



47a3|iff| + V^((i^f)'+7^3 



[KfY^ill, 
4 ll{{Kff+ll,) 

^^lal{\Kf\-^^a3f 

and the resonator mode amplitude is 

27a 



{b 



Vic 



{Ea), = 



Vi{\Kf\-V?.la3) 



Fluctuation 



(51) 
(52) 

(53) 



Fluctuation around the mean field solution are gov- 
erned by 



_d 
di 



(Ca\ 

'-a 
Cb 



W 



fca\ 

4 

Cb 

\4J 



fFa\ 

Fb 



(54) 



where the matrix W is given by 



W = 



/ Wi W2 W3 Wa 

Wi M^i* Wi W^ 

W5 We W7 Ws 

VWg* W^ Wi W^ 



(55) 



The mean field solution is assumed to be locally stable, 
that is, it is assume that all eigenvalues of W have a 
positive real part. 

We calculate below the statistical properties of the 
noise operators Fa and Fb- Let a(a;) be an annihila- 
tion operator for an incoming bath mode. In thermal 
equilibrium the following holds 

(a(c^))=0, (56) 

(a'' {uj)a{uj'))^n^S{uj-uj') , (57) 

(a {uj') flt (w)) = (ji„ + 1) ,5 (w - w') , (58) 

(a (w) a (w')) = , (59) 



where where 



1 .„.x w _ (Wi-ioj W2 



---^mi:-l^ (60) ^'^'^-(^ W| wt-^^)^ (^1) 

/3 = Xjk^T, fee is Boltzmann's constant and T is the VF^^ = { \ * ) , (72) 

absolute temperature. Usmg these expressions together V 4 "3 / 

with Eqs. (Uni, (HH), dMl), (EZD, (EHD, dMD and dSS]) yields _ /^ W^5 M^6 \ .7o^ 

the following relations *"" ~ i^ Wg* M^g* j ' ^ -^ 

(F, (..)) = (Fl (^)) ^ (F, (c.)) ^ (F,t (^)) =. , (61) W,, = [^ ^'^:'' ^Y~ ^UJ ) ' ^^^^ 

Multiplying the first equation by W~^ and the second 

{Fa (o.) Fa (o.')) - {Fl (^) Fi (^')) one by W-' leads to 

= (F, (a;) F, (..')) = (^/(^)^,^(^')) = 0, / CaH Vw-iw /" ^"H ^ w-^( ^'^H ^ 



(62) 



(F. (c.) Ft (.')) = 2r.^ (c. - c.') n.„ , (63) ( ^f^^"^^ ) +1^,-^W^,. ( /^"^^^^j) ) "'" ^ ^' ^"^ 

(F] (w) F, (w')) = 2ra(5 {uj ~ w') (n^„ + 1) . (64) or 



^- i Ft (-.) ; ' 



(75) 

^) ^ 

-). 
(76) 



(Fb (w) F,^ (cj')) = 276'^ (^ - w') '^-^ , (65) 



(77) 



and 

(F,t (c) F, (..')) = 27fa<5 (a; - u;') (n., + 1) , (66) (^ ^ (-a.) j " ^"^ (, F] (-c.) j + ^^^ (^ pi {-u) J 



where 



(78) 
where 
^'^ = ^^ + ^^^^''- ^ ^ (67) Xaa=(W^..-l^a.W^,^^l^..)-\ (79a) 

Is is important to note that the linearization approach / ti^-Itj^ \~i 



is valid only when the fluctuations around the mean field 



Xab = [Wba - WhbW-,'Waa) , (79b) 



solution are small. Unavoidably, however, very close to ^^^^ = (Wab — WaaWr^Wbb)~^ , (79c) 

the region where the system becomes unstable the flue- i -i 

tuations become appreciable, and consequently the lin- Xfcb — {'^bb — WbaWaa Wab) ■ (79d) 

earization approximation breaks down. ^j^^ ^^^^^^^ matrices W^-^ and W{;,^ can be expressed 



'bb 

as 



W^ Wi - iu 



{Wi - iuj) {W^ - iuj) - 


- \W2\^ 


f WI - iuj -W2 ^ 


1 


\ -W^ Wi-iio J 


1 



D. Transforming into Fourier space | '^i *^ '^2 

Waa 

In general, the Fourier transform of a time dependent 
operator O [t) is denoted as O (w) 

CXD 

O (t) = ^= f dujO (w) e-"^* . (68) ^ (Aai - Jw) {Xa2 - iw) 

V 27r J 

—00 

Applying the Fourier transform to Eq. (|54p yields / w* _ irj — M/ 

^-Ui(--)J+^"U:(-)J = Ui(-)J ' 

(69) 

^- ( 4 (-^i) ) - ^^< cU-^i) ) ^ ( fh'-I) ) ' ^ (^--)(^--) ' 

(70) 



(80) 



(VK7 - iw) (lyf - iuj) - 


- iWsl' 


f W^ - iuj -Wa ^ 




1^ -W^ Wt^iloj 





where we have introduced the eigenvalues 

AalAa2 = |T^l|'-|VK2|' , 



and 



Abi + Afc2 = 14^7 + W; , 
A6iAfc2 = IVK7I' - \Ws 



Using these relations one finds that 



(82a) 


/ 1 


\ 




(82b) 


Xaa- y 

/ 
1 V" 


1 




-Ea(Abi— Ab2) 

(Aai-iw)^ 

(B:)^(A61-A62) 


Bl{\ti-\b2) \ 




(Aai-iw)(A„2-ii^) 

Sa(A61-Ai,2) 


(83a) 


(A„i-iLj)(A„2-i'»j) 


(A„2-itu)^ / 


(83b) 




(Afei - 2w) (A 


,2 - iw) 



(91) 



E. Omega-Symmetric Matrix 

Let W (lj) be a 2X2 matrix, which depends on the real 
parameter uj. The matrix W (w) is said to be omega- 
symmetric if it can be written as 



W{uj) 



a (cj) b (uj) 
b* (— cj) a* (— cj) 



(84) 



where a (w) and b (w) are arbitrary smooth functions of 
Lj. Is is straightforward to show that if W is omega- 
symmetric then W~^ , W^ (transpose of W) and W'^ are 
all omega-symmetric as well. Moreover, if Wi and W2 
are both omega-symmetric then W1W2 is also omega- 
symmetric. Thus, it is easy to show that the suscep- 
tibility matrixes Xaa, Xab, Xba and Xbb are all omega- 
symmetric. 



F. The case where fl is small and Ka = 7a3 = 

To lowest order in Q one has 

Xaa - (1 - W-,'Wa,W,-,'W,ay' W, 

^ (1 + W-,'WabW,-,'Wba) W-,' 



-1 

aa 



Xab ^ 'W-,'WabW,-,' 



bb ' 



and 



Xba ^ -W,-,'WbaW-,' , 



Xbb = (1 - W-,'WtaW-,'Wab) ' W,-' 



(85) 



(86) 



(87) 



(1 + W,-,'WbaW-,'Wab) W^ 



-1 
bb 



Taking Ka = 7a3 = one has 



" aa \ Q 1 



Aa2 — iCJ 



Similarly W^^^^ can be expressed as 







w-l ^ Abi-iw 



bb 



Xb-2—i^ 



(88) 



(89) 



(90) 



i-B„ 



iS„ 



Xab — 



= -n 



(Aal— icj)(Abi— iw) (Aal— iw)(Ab2— itij) 

i-B* i-B* 



(Aa2-Jtj)(A(,i— iw) (Aa2-itj)(Ai,2 — itj) 



i-B;; 



Xba 



and 



Xbb 



= -0 



(Aal— icj)(A(,i— iw) (Aa2— i'»j)(Ai,i— icj) 

i-B* iSa 



(Aal-iw)(A(,2-Jtj) (Aa2— i'»j)(Ai,2-Jw) 



Abi —iuj 





i^2£;„ 



Ab2— itiJ 

(Aal-Aa2) 

(Abi-icj)^ 

(Aal-Aa2) 

(A(,i-iLj)(Ab2-ii<j) 



(Aal-Aa2) 

(A(,i-iw)(Ab2-i") 

(Aal-Aa2) 

(Ab2-ii<j)^ 



(Aal - »w) (Aa2 - iw) 



(92) 



(93) 



(94) 



To determine the stability of the mean field solutions 
the eigenvalues of W are calculated below for the present 
case to lowest nonvanishing order in fl. The matrix W 
can be expressed as 



W 



Aal 

Aa2 

Afti 

Ab2 



nv 



(95) 



where 



/ 





iBa 


iBa 








-ib: 


-iB 


^K 


tBa 








\-^B* 


-iBa 









V 



The two eigenvalues of interest for what follows are A;,i 
and Ab2, which approach the values Abi and A62 respec- 
tively in the limit fi — >■ 0. These eigenvalues are calcu- 
lated up to second order in J7 using perturbation theory 
(note that W is not necessarily Hermitian) 



Afei — Afci + Vt Ea 



\2 — Ah2 + ^ Ea 



1 



1 



Afel — Aal Afei — Aa2 

1 1 



Ah2 — Aal Ab2 — AaS 



, (96a) 
(96b) 



Thus by using the relations 



the notation 



Afei = Xl-^ = iujb + lb , 



LUb 
la 

9^ — , 

CJb 



(97a) 
(97b) 



(98a) 
(98b) 



and by assuming also that ib ^ ^b one finds that 



Xbi = iw& 1 



+ 7b 1 + 



2n''Ea 


2d(l- 


Q 


^-5^) 


2 
^6 


[d + if + .g2 




(rf-l)Vg2 


2Q?Ea 


4dg2 ^ 



lalb 



{l + df 



{l-df + g^ 



(99) 



and Ah2 = Aj;^. 

For the present case {Ka ~ lai ~ 0) one finds using 
Eqs. dMI) and ^ that {Ea)^ (the value of £'q at the 
onset of bistability) is given by 



{Ea), 



lal^b 



(100) 



V302 ■ 
In terms of (Ea)^ the real part of Xbi can be expressed 



as 



Re 



M.i 



76 



V3(Sa),7^ 



T(d,ff), (101) 



where the function T {d, g), which is plotted in Fig. [TJ is 
given by 



r{d,g) = 



Mg^ 



{l + df+g-^ {l-df+g'^ 



4g^d 



4<?2 + (g2 _ m_ ^2y 



(102) 



For any given value of g the function T obtains a maxima 
at d = do and a minima at d = —do, where 



1 



do^-yS-Sg^ + G^g^+g^ + l. (103) 

The mean field solution is stable provided that 
Re ( Afci J > 0. Hopf bifurcation occurs when Re ( Xbi ) 
vanishes. 



0.5 



-0.5 




FIG. 1: The function T (d,g). 



IV. INTEGRATED SPECTRAL DENSITY 

In general consider an operator c (w) that can be ex- 
pressed in terms of a noise operator F (w) and a suscep- 
tibility matrix xi'^) as [similarly to Eqs. (iTTl) and ((78)) ] 









(104) 



where F (w) satisfy [similarly to Eqs. dH]), dM]), (pg]) . 
del, dM]) and dMl)] 



(105) 



(F(c.)) = (Ft(c.))=0, 
(i^ (w) F (cj')> = (F^ (oj) F^ (w')) = 
{F{Lj)FHij'))=2TS{ij-uj')n^^ , 



(106) 



(107) 



and 



by 



(Ft (^) i^ (^')) ^ 2T5 {uj - uj') (n^„ + 1) . (108) 
The homodyne detection observable X (w) is defined 



X (uj) = e'*'^°c (cj) + e-"^'^°c^ (-w) 



(109) 



The frequency auto-correlation function of X is related 
to the spectral density Px (w) by 



{X^ {lu') X (uj)) = Px H Siuj- uj') 



(110) 



Assuming that xi'-^) is omega-symmetric, it can be ex- 
pressed as 



X(w) = 



a (uj) b [oj) 
b* i-uj) a* (-w) ; ' 



(111) 



8 



where a (w) and b (w) are arbitrary functions of uj. Thus, 
by calculating the term (^X^ (w') X (w)) one finds that 



2P = M+ (w) coth — ^ + A/_ (w) 



(112) 



where 



\a{^^)\' + \b{uj)\' + \a{uj)\' + \b{-Lo)\' 



+ Re [e^'*^" (a (-w) fe (w) + a (cj) & (-w))] 



M+ (w) 

and 

M_ (uj) 

+ Re [e^*"^^" (-a (-w) 6 (w) + a (w) 6 (-cj))] . 

(114) 
The integrated spectral density (ISD) is thus given by 

oo 



(113) 



|a(-^)|^-|6(c.)|V|a(^)|V|6(-c.)|' 



the case where resonator b has high quality factor is as- 
sume. For this case, which is experimentally common, 
the following is assumed to hold 7^ ^ Wf, and 7fc ^ 7a- 
As can be seen from Eq. (|118l) . for finite driven ampli- 
tude Ea the ISD of Xii can deviate from the equilibrium 
value of coth {j3fvjJi,/2). 



VI. DECOHERENCE 

The Hamiltonian of the system ([2|) is formally a func- 
tion of Ab and A^, that is % = "H (A;,, AH. Consider 

resonator 6 in a superposition of two coherent states |ai) 
and \a2)- In therm of the operator V, which is given by 

V = 'H{a2,a*2)-n{ai,al) , (119) 

the decoherence rate 1/r^ can be expressed as [43| 



1 



/ Auj Px{oj) = 2rV"coth- 



1 



dw (V(O)V(w)) , 



(120) 



(115) 



where 

c 

y = / dw M+ (w) 



where 

V (i) = V (<) - (V {t)) , (121) 

and where V (w) is the Fourier transform of V (i) 
1 



V{t) 



/2ti 



dw V (w) e 



(122) 



duj 



\a (w)l' + |fe (-cj)l' + 2Re (e^^^^^a (w) fe (-w)) 



(116) 



ISD OF Xft 



We calculate below the ISD of the homodyne observ- 
able Xf, (w), which is given by 

Xb (w) = e"^^^Cb {u) + e-">'^°cl {-uo) , (117) 

for the case where fi is small and Ka = 7a3 = 0. As 
can be seen from Eq. (|75|) . it has two contributions due 
to the two uncorrelated noise terms Fb (w) and Fa (w). 
The calculation of both contributions according to Eq. 
(jllSp is involved with evaluation of some integrals, which 
can be performed using the residue theorem. To further 
simplify the final result, which is given by 



' Using Eqs. (|35p and ([77| together with the notation 
<^a = "2 - "1 = |<5a| e' , (123) 

one finds to lowest order that 

^ = UaFa {u:) + U:fI i-uj) + UbFb ic,) + U:F^ {-cu) , 

(124) 
where 

Ua = 2n cos e{B:{Xaa) II +BaiXaa)2l) , (125a) 

Ub = 2n COS 6(3: {xab)n + Ba iXab)2i) + ^e-'" . 

(125b) 

Furthermore, with the help of Eqs. ([55)) . (jM)) . (|55|) and 
(pS)) the decoherence rate becomes 



oo 






7a 7& 

»'^a 2g^ 

7a7ft (l-d)^-t-g2 



■ coth 



2\5^f (Tapaf coth ^+jb\Ub\' coth ^ 



(126) 

Note that for fi = the decoherence rate reproduces the 
value given by Eq. ([T]). 

For the case where H. is small and Ka = jas = one 
finds using Eqs. ^ and ^ that 



(118) 



\Ua\ 



An^Ea cos2 e 



(127) 



and 



\u,\' = i- 



where 



-iVL^Ea [cos a + cos {29 + a)] d 

<P +g2 



^hn 



ii) 



tan 



-1 7^ 

UJb 



(128) 



(129) 



In what follows we restrict the discussion to the case 
where 9 = 0, for which the two coherent states |ai) and 
1 02) have the same momentum. For this case, which 
is the assumed case in some of the pubhshed proposals 
for observation of quantum superposition in mechanical 
systems d, [3, |4y], up to first order in ^t/^b one has 



\Ub\' 



1 



in^E, 



(p 



(130) 



Using these results together with Eq. (|126p one finds 
that 



- = 27fc|<5„|' 



4172^;, 



(P 



coth 



fifu^:b 



la ^^^Ea 1 



7b ^E 



(p 



■ coth 



PhuJa 



(131) 



The first term in Eq. (|13ip represents the contribution 
of the thermal bath that is directly coupled to resonator 
h to the dephasing rate. This contribution can be either 
enhanced {d > 0) or suppressed {d < 0) due to back- 
reaction effects. On the other hand, the last term in Eq. 
(|13ip [compare with Eq. (71) of Ref. ji^] represents 
the direct contribution of the driven resonator a. This 
contribution can be understood in terms of the shift in 
the effective resonance frequency of resonator a between 
the two values corresponding to the two coherent states 
|ai) and [02) (see Ref. J47|). 



— = 2-fb\Sa\ -^ , 



(133) 



where 



IiSD ^ ^ _ n-'EaT{d,g) f _ aJbOg {1 + df + gA 

T lalb y i^a id J ' 

(134) 

Tu ^ , 4n^Ea d (^ , LOblaQa l\ .^„,. 

^ = 1 + ^.2 ^2 , „2 1 1 + — 7T— J ) ' (135) 



d^ + 92 



Wa7fc d 



and 



Ha = coth 

2 2 



(136) 



In terms of (Ea)^, which is given by Eq. (|100[) . one thus 
has 



T 



and 



= 1- 



Eg ujbT{d,g) 

(Ea)^ VSjb 



1- 



iObOa {I + dy + g^ 



3d 



Td 
T 



4.E„ 



dg 



VS{Ea)J^+9^ 



Wa76 d 



(137) 



(138) 



These results are valid only to lowest order in fJ, how- 
ever they may be used in some cases to roughly estimate 
the lowest possible values of Tisd and Tq. As can be 
seen from Eqs. (|137p and (|138p . the effective tempera- 
tures TisD and Td rnay take considerably different values. 
This fact should not be considered as surprising since the 
system is far from thermal equilibrium and since the un- 
derlying mechanisms responsible for ISD reduction and 
for suppression of decoherence are entirely different. In 
what follows, we choose the parameters d and g such that 
the largest reduction in effective temperature is achieved 
for a given Eg , and use these values to estimate the lowest 
possible effective temperatures. 



VII. DISCUSSION 



A. Optimum ISD Reduction 



We have considered above the case where fl is small. 
Kg = 7a3 — and jb ^ Wb. In addition, we have as- 
sumed that 7a <C 76 in order to obtain the ISD of Xb, 
which is given by Eq. (|118|) . and we have assumed the 
case 9 = to obtain the dephasing rate, which is given by 
Eq. ()13ip . Furthermore, consider for simplicity the case 
of high temperature where phub <C 1. For this case Eqs. 
piSp and (|13ip can be written in terms of the effective 
temperatures Tjsd and Td 



00 
^/d.Px,(.) = 



2fcBriSD 
hu!b 



(132) 



For the case of ISD reduction, we consider the case 
where the term that is proportional to Qg in Eq. (I137p . 
namely the term which represents the contribution of the 
thermal baths that are directly coupled to resonator a, 
is relatively small, namely the case where uJbOg <^ ujg. 
This condition is expected to be fulfilled for the typical 
experimental situation. Most efficient ISD reduction is 
achieved by choosing the parameters 17 <C 1 and d = 1, 
for which the term T (d, g) obtains its maximum possible 
value T = 1 (see Fig. [T]). For this case Eq. (|137p becomes 



TlSD 



^b Eg 

V37b (Eg), 



4:UJbQg 



(139) 



10 



By taking 



En - (^a)c = (-E^a)lSD ' 



LOb 



(140) 



Eq. p39p yields the lowest possible value of Tisd , which 
is denoted as (TisD)min 



(^isd)i; 
T 



4^fc9a 

3wq 



(141) 



As was mentioned above, the above discussion is based 
on the approximated result Eq. (I137p . which expresses 
TisT) to lowest nonvanishing order in f2. Such an expan- 
sion apparently suggests that the noise contribution due 
to the thermal bath that is directly coupled to resonator b 
can be altogether eliminated, leaving thus only the noise 
contribution of the thermal baths that are directly cou- 
pled to resonator a as a lower bound imposed upon Tisd 
[see Eq. (I14ip ]. Obviously, however, higher orders in Q, 
have to be taken into account in order estimate (TisD)i„in 
more accurately, as was done in Ref. [4^ , where Tisd was 
expanded up to forth order in fi. 



B. Decoherence Suppression 

For the case of decoherence suppression, on the other 
hand, the term that is proportional to Qa in Eq. (|138p 
is not necessarily small for the common experimental sit- 
uation. We therefore chose the optimum values of the 
parameters d and g for the more general case. Using the 
notation 



D 



OJa^b 



Eq. (IT^ reads 






where 



fid,g,D) 



dg 



D 

g- \ d 



(142) 



(143) 



(144) 



In general, the minimum value of the function / (d, g, D) 
for a given g > and a given Z? > is obtained at 



d,„ = -D- VD^+g^ , 

and the minimum value is given by 

1 /tan-i j. 
f{d^,g,D)^~-t8ini ^ 



(145) 



(146) 



The lowest value of f {djn,g,D) is thus obtained in the 
limit D ^ g, for which one finds that d^ = —g and 



f {dm,9,D) — —1/2. Therefore, one concludes that the 
largest reduction in Tq for a given Ea is obtained when 



Wa7h 



<1 



(147) 



and when d — —g. For this case Eq. (I138P becomes 



lb 
T 



1 



Er, 



ViiEa) 



(148) 



This results indicates that even when all parameters are 
optimally chosen such that the largest reduction in Td is 
obtained for a given Ea, no significant reduction in Td 
is possible unless Ea becomes comparable with (Ea)^. 
Note, however, that in our analysis of the present case 
the effect of nonlinear bistability has been disregarded. 
This approximation can be justified for the case of ISD 
reduction since, as can be seen from Eq. (|140p . opti- 
mum reduction of the ISD can be achieved well below 
the bistability threshold provided that 7^ ^ Wf,. On the 
other hand, Eq. (|148p indicates that optimum suppres- 
sion of decoherence can be achieved only very close to 
the bistability threshold. In this region, however, our 
approximated treatment breaks down and Eq. (|138p be- 
comes inaccurate. 

To calculate Td near the bistability threshold we thus 
numerically evaluate the dephasing rate given by Eq. 
(|126p without assuming that fl is small or Ka = 7a3 = 0. 
As before, we take 9 = and consider for simplicity the 
case where (3fiuJb ^ 1, for which the effective temperature 
Td is given by 



^ = \Ub\' + ^^^^^^^^^^ea 



(149) 



Figure [5] shows an example calculation of the parame- 
ters \Ub\ and \Ua\ and the ratio Td/T near bistability 
threshold of the system. The ratio Td/T is shown for the 
case where fShtOa ^ 1- The set of system's parameters 
chosen for this example is listed in the caption of Fig. 
[21 The stability of the mean filed solution is checked by 
evaluating the eigenvalues of the matrix W. The dot- 
ted sections of the curve Td/T indicate the regions in 
which the solution is unstable (where at least one of the 
eigenvalues of W has a negative real part). Near the 
onset of bistability point [see panel (c4) of Fig. [5] and 
near jump points in the region of bistability [see panel 
(d4) of Fig. [2] the ratio Td/T may become relatively 
small. This behavior can be attributed to critical slow- 
ing down, which occurs near these instability points [43] . 
On the other hand, in the vicinity of these points the 
solution becomes unstable [see the dotted sections of the 
curve Td/T in panels (c4) and (d4) of Fig. [2]. When 
the unstable region is excluded one finds that no signifi- 
cant reduction in the ratio Td/T can be achieved for this 
particular example (the lowest value is about 0.5). 

In the previous example the mean-field solutions be- 
come unstable close to the onset of bistability. This be- 
havior prevents any significant suppression of decoher- 
ence, namely, the ratio Td/T could not be made much 
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FIG. 2: The factors \Utf and \Ua\'^ and the ratio Td/T. The 
driving amplitudes in columns a, b, c and d are bp/ (bp)^ ~ 
0.01, 0.8, 1 and 1.3 respectively. Other system parameters 



are 0,/uja 



10" 



i^b/oJa 



10" 



'yb/ojb 



10-^ Ka 



2 X 2n^LJi/ (ujl+'yi), 6» = 0, -Jal/uJi = 10^ 7a2/7al = 10"^ 

and 7„3 = 0.1 x \Kf^\ /Vs. The ratio Td/T, which is plotted 
in the forth row, is shown for the case phua ^ 1. The dotted 
sections indicate instability. 



smaller than unity. To overcome this limitation the pa- 
rameter (Ea)^, which is given by Eq. (j53p . has to be in- 
creased without, however, increasing the coupling param- 
eter fl. We point out below two possibilities to achieve 
this. In the first one, the parameter Ka is chosen such 
that Ka — 2n'^uji,/ (ujI + 7 ^) , and consequently K^^ be- 
comes very small [see Eq. (|45p ]. In the second one, which 
is demonstrated in Fig. [3] below, the nonlinear damp- 
ing rate 7^3 is chosen very close to the largest possible 
value of \Ka^\ /\/3 for which bistability is accessible [see 
inequality ((50|)]. As can be see from Eq. ([53|). both 
possibilities allow significantly increasing the parameter 
(Ea)^. For the example shown in Fig. [3] below, the value 
7a3 = 0.99 \Ka^\ /-\/3 is chosen and all other parameters 
are the same as in the previous example (see caption of 
Fig. [2]). As can be seen from panels (c4) and (d4) of Fig. 
[21 much lower values of the ratio Td/T are achievable 
in the present example (a lowest value of about 0.02 is 
obtained at the edge of the region where the solution is 
stable). This improvement can be attributed to the sta- 
bilization effect of the nonlinear damping. It is impor- 
tant to point out, however, that implementation of any 
of the two above mentioned possibilities require that the 
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FIG. 3: The factors \Ubf and \Uaf and the ratio Td/T. In 
this example 7a3 = 0.99 x \Ka^\ /\/3 whereas all other pa- 
rameters are the same as in the previous example [see caption 
of Fig. (0]. 



nonlinear parameters of resonator a (Ka and/or 7^3) can 
be accurately tuned to the desired values. Such tuning 
of nonlinear parameters can possibly becomes achievable 
by exploiting effects arising from thermo-optomechanical 
coupling j49| . 
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